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ABSTRACT 

A method is presented for the expeditious analysis of performance 
of a solution gas drive reservoir. This rapid analysis is accomplished by 
the repeated solution of the Muskat equation on the IBM- 701 computer. In 
each solution, one or more reservoir parameters are varied to determine 
the effect on the problem results. A small California reservoir is used to 
demonstrate the procedure. Graphs, illustrating the results for each 
parameter range, are presented. 

For the sample reservoir only production, depth-pressure, and 
electric log data were available. Average values of PVT, viscosity, and 
permeability were developed from published correlation in California Oils, 
and entered into the Muskat prediction computer program. ^ Fifty predictions 
were run, and, of these, twenty-three solutions are graphically illustrated. 
From the graphs, a set of parameters was selected which best define the 
sample reservoir. 

The results of this analysis indicate that more than one prediction 
is necessary to establish the reservoir parameters. By continuously ob- 
serving the predictions, the reservoir engineer is able to decide when a 
reasonable approximation of the reservoir's parameters has been made. 

He must know, from experience and his knowledge of reservoirs, which 
set of parameters best defines the reservoir. 
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INTRODUCTION 

For the last decade, Muskat‘s equation^ for the depletion history 
of a solution gas -drive reservoir has been a reliable tool of the reservoir 
engineer. It undoubtedly would be used more extensively, but for the vast . 
amount of time and labor required to complete one Muskat reservoir pre- 
diction. (An estimated sixteen man hours per prediction). Often, some of 
the data which the reservoir engineer must use are doubtful. It is desirable, 
therefore, that the Muskat prediction be run a number of times, using an 
array of values for the doubtful data. However, the man hours required for 
such a procedure would soon become prohibitive. The obvious solution to 
this dilemma is the high speed digital computer. Muskat's material balance 
equation, expressed in differential form, is ideally suited to digital computer 
solution. 

The purpose of this paper is to illustrate how a series of Muskat pre- 
dictions, obtained in a relatively short time, provided the reservoir engineer 
with a production analysis of his reservoir. 

A small California reservoir was used to demonstrate this procedure. 
Data for the reservoir were provided by the Standard Oil Company of California. 
These data included electric logs, production records, and a limited amount 
of depth pressure data. No PVT data or core data were available. Porosity 
and water saturation data for a similar zone, in a nearby reservoir, were 
used to make initial estimates. 
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THEORY 



The IBM- 701 computer solution to the Muskat Prediction, once pro- 
grammed, may be run, with slight alteration to the variables, until the pre- 
diction data and the production data coincide. When this occurs, it may be 
assumed that the prediction, extended beyond the current produced life of 
the reservoir, will be an accurate estimate of the reservoir's future, pro- 
vided the operating procedure of the reservoir is unchanged. The computer 
makes a single prediction in about three minutes. However, if three pre- 
dictions are run in tandem, with the same imput quantity varied in each (e. g. 
a different size gas cap in each run), the three predictions will be completed 
in less than seven minutes. 



The IBM-701 program herein utilized was developed by J. E. Warren 

O 

and T. D. Mueller' 1 of the California Research Corporation at La Habra, 
California. The program includes the equations listed below: 



For production above the bubble point pressure: 
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where: 

Np » Cumulative oil production in bbl of STO/acre ft. 

Gp a Cumulative gas production in MCF/acre ft. 

R * Produced gas -oil ratio in SCF/STbbl. 

All other symbols are defined in the appendix, page i. 

The computer program is designed to solve these equations, and 
to record their solution as illustrated by Figure 1. 

To adapt the computer program to a particular reservoir, the PVT 
data, viscosity data, and permeability data for that reservoir are introduced 
into the program. To consider a new reservoir, one simply substitutes new 
data. 

RESERVOIR HISTORY 

From the electric logs provided, a sand count was made of the pro- 
ducing sand, and a contour map was developed, (see Figure 2). It was 
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apparent, from the production data and from numerous cross sections of 
the producing zone and the two zones above it, that the reservoir provided 
was actually three reservoirs, separated from one another by fault zones. 

It was decided that only one of these reservoirs. Reservoir B, (illustrated 
in Figure 2) would be used to demonstrate the IBM- 701 prediction procedure. 

The contour map of Reservoir B shows that the reservoir contains 
nine wells. Of the nine wells drilled, only two are currently producing. 

Well #5 was the first well brought in (November 1949) with an initial pro- 
duction of 400 bbls/D. It produced until December 1950 (a total of 82, 000 
bbls), when it was shut in with a maximum gas-oil ratio of 7500. Well #18 
was completed in January 1950, initial production 9 bbls/D. It produced 
until December, 1951 when its maximum gas-oil ratio was 6700. It produced 
a total of 2, 000 bbls of oil. Well #1 was drilled in June 1950, and came in a 
gas well. It was immediately shut-in. Well #8 was completed in November 
1950; its initial production was 380 bbls/D. It is currently producing 31 bbls/D. 
Its gas-oil ratio has climbed from 280 to 4, 000. It has produced a total of 
185, 000 bbls. of oil. Well #9 was drilled in January 1953. It came in a gas 
well and was shut in after producing 3, 600 MCF of gas and 33 bbls. of oil. 

Well #12 was completed in January, 1956 with an initial production of 49 bbls/D. 
It ha6 produced 22, 800 bbls. to date and is currently producing 51 bbls/D. This 
well initially produced 118 bbls/D of water, and is now producing 272 bbls/D. 

It is the only well in the field that has produced a sizable amount of water. 

Wells #2, #4, and #14 were dry holes. As illustrated in Figure 2, the entire 
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northern part of Reservoir B, which includes wells #5, #18, #1, #2, and #9 
R/D, is a gas cap. There is also evidence of water influx in the eastern 
edge of the reservoir near well #12. 

The total oil produced to date from Reservoir B (July 1957) is 507, 356 
bbls. of stock tank oil. This oil has an average gravity of 35° API. 673, 457 
MCF of gas has been produced. 

PREDICTION PROCEDURE 

Production data provided an average API gravity of 35° for stock 
tank oil, an average gas gravity of 0. 73 for the reservoir gas at standard 
stock tank conditions, and an average reservoir temperature of 135° F. 

o 4 

Standing's correlation charts were used to obtain PVT data, and Beal's 

correlations were used for the viscosity data. PVT data and viscosity data, 

as a function of pressure, are included as Graph A. 

From the sand count of Reservoir B, an isopachous map was constructed 

(Figure 3). A volumetric estimate of sand in place was computed to be 2,520 

acre feet. Porosity and connate water saturation were estimated to be 0. 35 

and 0. 45 respectively. These were the average values of twenty-two porosities 

and water saturations taken from sand in an adjacent reservoir. The adjacent 

sand was considered to be similar. Using the equation: 

N a (7756) (acre ft. ) U ) (1 - S y) , 

B oi 

an initial estimate of oil in place was computed to be 2. 88 million stock tank 
barrels. (Sample calculation, page ii in appendix illustrates this procedure). 
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GRAPH A 
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This estimate was substantiated by the Schilthuis Material Balance Equation 
when an allowance was made for gas production from the gas cap. A second 
calculation of N was made using new values of porosity and connate water 
saturation, <t> = 0. 25 and S w =0. 30. These estimates were taken as average 
for California sandstone. The initial estimates are on the extreme end of 
the range for California sandstones. The new values of porosity and water 
saturation give a value of N equal to 3.14 million barrels of stock tank oil. 

The change in N is comparatively slight. 

A curve of oil saturation vs. relative permeability, (Graph B), was 
developed by solving the equations listed below: 

7. S 0 =(l-S w )(l-^)f^ 

8. k 9/k 0 = (R-R S ) M9 /m 0 7g 0 

Static bottom hole pressure, cumulative oil and gas production, and 
instantaneous gas-oil ratios were plotted vs. time in years, (see Graphs C 
& D). From the correlations established, R and AN (bbls. of STO) were 
taken as functions of pressure. Dividing A N by N gave A N/N as a function 
of pressure. Since two N's were computed, two cases. Case I, where <t> - 
0. 35 and S w x 0. 45 and Case II, where <b = 0. 25 and S w c 0. 30 were developed. 
PVT and viscosity data, as functions of pressure, were taken from Graph A. 
Equations (7) and (8) were solved as functions of pressure in decrements of 
100 psi from 1500 to 600 psia. From the correlation of relative permeability 
and oil saturation for the same decrement of pressure, two curves of Kg/K Q 
vs. S Q for Case I and Case II. Both curves were extended into the region of 
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OIL SATURATION, PERCENT 



GRAPH B 
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higher gas saturation by continuing them parallel to the curve developed 
by Arps^ for the average value of relative permeability in sandstones. Once 
the PVT, viscosity, and relative permeability data were compiled, they were 
included in the computer program. Experience with the computer program 
indicates that most functional data, used by the computer, must be mathema- 
tically smoothed. The smoothing of the PVT and viscosity data was accom- 
plished by third order polynomial equations. A seventh order equation is 
usually required for Kg/K Q functions. However, in this particular case, 
the curve from Arp's average was sufficiently smooth to enter the program 
directly. 

PREDICTION RESULTS 

It was anticipated that the predicted results could be brought to 
superimpose upon the observed production data. This was to be accomplished 
by a system of altering the problem variable until the desired results were 
reached. 

For the first prediction, it was decided to investigate the Case I 
parameters. The following initial conditions were set: 

P^ = 1500 psia, bubble point pressure = 1500 psia 
S w = 0. 45, <i> = 0. 35, m * 0. 1, S 0 i* - 0. 15 
The predicted results are illustrated in Figures 4. 0 and 4. 1. * 



^Predicted data was illustrated in all figures by a solid curve. Observed 
data was shown as circled points. 
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Figures 5. 0 and 5.1 show that lowering the bubble point causes the 
predicted data to begin to approach the observed data. The best approxirna - 
tion of % A N/N is for a BPP* 900 psia, but the declining slope after P = 

900 is too gentle. 

After this prediction was made, the effect of varying porosity, water 
saturation, viscosity, gas cap and initial pressure was investigated. None 
of these variations produced a unique solution. Table 1. 0 lists the parameters 



investigated. 
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Very little agreement is shown between observed and predicted 

data. 

To investigate the effect of bubble point on the prediction a second 
run was made. This was a tandem run with a different bubble point used 
for each prediction. The following initial conditions were set: 

Pj = 1500 psia.BPP s 1500, 1200, 900 psia, S w = 0.45, <j> = 0.35, 
m = 0.1, S Q £ * 0. 15 



The predicted results are illustrated in Figures 5. 0 and 5.1. 



TABLE 1. 0 
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The Case II parameters were next investigated. For the third pre- 
diction, the following initial conditions were set: 

* 1500 psia, BPP * 1500 psia, S w * 0.30, $ * 0.25, m =0.5, 

S .* = 0.15 
oi 

The predicted results are illustrated in Figures 6.0 & 6. 1. 
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Here the GOR curve showed slight agreement. No agreement was 
apparent for oil production. 

The fourth prediction was used to investigate the effect of gas cap 
size. The following initial conditions were set: 

Pi * 1500, BPP * 1500 psia, S w * 0.30, 4 * 0. 25, m = 0.1, 0.5 
& 0.9, S Qi * = 0.15 

The predicted results are illustrated in Figures 7.0 and 7.1. 
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The fourth prediction indicated that the gas cap size does have a 
significant effect upon gas and oil produced. It suggests that m = 0. 1 is 
probably the more likely size for this reservoir. There was no agreement 
between predicted and observed data for oil production, but predicted and 
observed gas-oil ratio compared well. 

The fifth prediction was used to demonstrate the effect of bubble 
point upon production. For this prediction the following initial conditions 
were set: 

Pi - 1500 psia, BPP = 1500, 1200, 900 psia, S w * 0.30, <t> =0.25 
m =0.1, S oi * = 0.15 

The predicted results are illustrated in Figures 8.0 and 8.1. 



18 . 




The fifth prediction showed that lowering the bubble point pressure 
caused the predicted oil curve to approach the curve of the actual oil produced. 
At a bubble point of 900 psia, the predicted oil curve superimposed upon the 
actual produced oil curve. However, the corresponding GOR curve atBPP * 
900 psia was too low. This may be explained by examining the history of 
Reservoir B. The excessively high gas oil ratios which occur at various 
times throughout the reservoir's history, strongly suggest that gas was 
produced from the gas cap. Since the Muskat equation does not account 
for gas production directly from the gas cap, the computer program was 
not designed to reflect it. It is concluded that the difference between GOR 
computed and actual GOR is due to the gas produced directly from the gas 
cap. The prediction indicates that slightly over half of the gas produced 
from the initial pressure of 1500 psia to 600 psia was produced directly 
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from the gas cap. A new estimate of actual GOR curve was made excluding 
the gas cap production. This new GOR will be reflected on the remaining 
curves. The fifth prediction suggests that a bubble point of approximately 
900 psia is most likely for Reservoir B. Using this new bubble point results 
in a new value of N, based on a change in B 0 i. The change is slight, however, 
and not sufficient to alter the Kg/K 0 vs. S Q curve. 

It should be noted on figure 8. 0 that the recoverable oil at abandonment 
pressure of 200 psia was roughly the same for all three bubble points. 

To substantiate the earlier decision to use a gas cap size m * 0.1, 
in the sixth prediction, a variable gas cap was once more employed; this 
time with the bubble point pressure at 900 psia. The following initial 
conditions were set: 

Pi = 1500 psia, BPP e 900 psia, Sw * 0.30 

<t> * 0. 25, m = 0.1, 0. 5, & 0. 9, S 0 j * 0. 15 
The predicted results are illustrated in Figures 9.0 and 9.1. 
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Once again, an m = 0. 1 appeared to be the best estimate for gas 
cap size in Reservoir B. 



To investigate the influence of viscosity, the seventh prediction was 
made with the following initial conditions: 

Pj = 1500 psia, BPP * 900 psia, S w * 0. 30, * 0. 25, 

m *0.1, S oi * * 0. 15, and 0. 8 M, 1. 0M & 1. 5M where M * 

The predicted results are illustrated in Figures 10. 0 and 10. 1. 
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As one would expect a higher viscosity results in less recoverable 
oil and a consequent higher GOR. The results indicated that the initial 
estimate of oil viscosity ( M * jij x 1. 0 ) is the best approximation for 
the actual production data. 

In the eighth prediction, the influence of porosity was investigated. 
Changing porosity changed the value of N, the initial oil in place. It did not 
change the value of % A N/N because porosity cancels out of this ratio. 
However, in the calculation of oil and gas produced per acre foot, porosity 
is reflected. Converting the acre feet production to predicted reservoir 
production permits a comparison with actual production, Graph D. 

The change in N for each value of porosity resulted in a very slight 
shift of Kg/K Q along the S Q axis (Graph B and Equation 7). The shift was 



22 . 



so slight that Arp's average curve was not varied. 

The eighth prediction initial conditions were set as follows: 

- 1500 psia, BPP * 900 psia, S w = 0.30, i> * 0.20, 0. 25 & 0. 30, 
m =0.1, S 0 j* * 0.15 

The predicted results are illustrated in Figures 11.0 and 11.1. 





This prediction is compared with actual production data illustrated on 
Graph D. The production data is not altered by dividing it by N. It is concluded 
that for Reservoir B, a porosity of 0. 25 is the only porosity which will give a 
prediction that satisfies actual production data. This condition holds so long 
as Kg/K 0 curve was shifted along the S Q axis for each change in S w . This shift 
was not due to change in N, but due entirely to the change in S w . 



The ninth prediction was made to investigate the effect of connate water 
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saturation. With the porosity held constant at 0. 25, the change in S w caused 
a slight change in N. To illustrate how this change in N caused a shift of 
Kg/K Q curve along the S Q axis, a sample calculation was made. N is expressed 
in millions of stock tank barrels. When positioning the Kg/ K q curve on the 
S Q axis, S Q was obtained by solving Equation (7) at various values of pressure. 
Example: when P * 700 psia. 

So ■ (1 - S w > (1 - > |p 

at N * 3. 12 when <t> = 0. 25 and S w ~ . 30 
S Q = (. 7) (1 - - 41 ) (. 865) = 0. 56 

<3 • 1 Ct 

if N = 2. 9 for <t> = 0. 25 and S w * 3.5 

S n * (. 65) (1 - • 41 ) (. 865) « 0. 515 , a shift of 4. 5%. 

2.9 

However, if only the change in N is reflected, 

N * 2. 9 for <t> = 0. 25 and S w * 0. 30 

S o = (• 7) (1 - ) (.865) = .555 

or the shift of the K ff /K Q curve along the S D axis for a change in N of 220, 000 

o 

bbls. of STO is only 0. 5%. It was concluded that the small changes in N, due 
to changes in porosity and water saturation, were not sufficient to shift the 
Kg/K Q curve, but the changes in S w did move the curve along the S Q axis a 
distance equal to that change. Consequently, it follows that Kg/K Q vs. S Q 
curve, taken as Arp's average for sandstones, is a fixed parameter except 
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when the value of S w is varied. 

Initial conditions for the variable connate water saturation prediction 
were set as follows: 

p. = 1500 psia, BPP = 900 psia, S w 3 0. 25, 0. 30, and 0. 35, <t> = 0. 25 
m * 0.1, S 0 i* * 0. 15. 

The predicted results are illustrated in Figures 12.0 and 12.1. 




Figures 12.0 and 12.1 illustrate the effect of the new Kg/K 0 vs. S Q 
curves on the prediction. These curves, which were extended by Arp's 
average, satisfied the observed data only when S w = 0. 30. 

The tenth prediction was made to investigate the effect of initial 
pressure. Initial conditions were set as follows: 

p. = 1500, 1400, 1300 psia, BPP 3 900 psia, S w = 0.30, <t> = 0.25, 
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m * 0.1, S 0 j* * 0. 15. 

The results were not of sufficient difference to merit plotting. The oil 
recovery at abandonment pressure of 200 psia for each Pj was: 



P^ * 1500 psia, % 


AN/N * 23.57 


1400 psia, % 


AN/N » 23.49 


1300 psia, % 


AN/N = 23. 39 



It was concluded that an initial pressure from 1500 to 1300 psia will 
have only slight effect on recoverable oil. 

The eleventh prediction illustrated the effect of changing S 0 i*. The 
following initial conditions were set: 

Pi * 1500 psia, BPP = 900 psia, S w = 0. 30, <t> * 0. 25 
m = 0.1, S Q i* = 0. 10, 0. 15, and 0. 20. 

The results indicated that, within the average limits of S Q i* for a 
California solution-gas drive reservoir, S Q i variation has little or no effect. 

It was concluded that the following parameters describe Reservoir B: 

Bubble point pressure 900 psia 

Initial reservoir pressure 1300 - 1500 psia 

Connate water saturation 0. 30 

Porosity 0. 25 

Ratio of initial reservoir free gas volume to initial reservoir oil 



volume 0.1 
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Oil saturation in the gas cap 0. 15 

PVT, viscosity, and relative permeability - - - as described in Graphs 
A and B, respectively 

Initial oil in place in Reservoir B was estimated to be 3.12 million 
barrels. Recoverable oil was estimated to be 730, 000 barrels of stock tank 
oil or 23. 57% of initial oil in place. 

DISCUSSION OF RESULTS 

The prediction results* have illustrated how quickly the reservoir 
engineer may bring predicted data to superimpose upon actual production 
data. By the fifth production, it was reasonably apparent that recoverable 
oil was about 23. 6% of oil in place. The remaining predictions simply re- 
affirmed the fifth prediction by illustrating that the Case II values, selected 
for the variable data, gave solutions which more nearly approximated the 
actual production data. Such will not always be the case. However, when 
the reservoir engineer runs the Muskat prediction in tandem, with three 
different values for some doubtful input, he will very rapidly pick up a 
trend which will enable him to close in on the parameters which specify his 
reservoir. 

Reservoir B was not the best possible reservoir to illustrate this 
problem. It is, however, typical of the situations which the reservoir 
engineer continuously encounters. 

From the history of Reservoir B, it was obvious that gas was often 



♦See Table 2. 0. 



TABLE 2. 0 



Summary of Parameters Used 



Prediction 

number 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Initial 




















1500 


Pressure 


1500 


1500 


1500 


1500 


1500 


1500 


1500 


1500 


1500 


1400 






















1300 


Bubble 




1500 






1500 












Point 


1500 


1200 


1500 


1500 


1200 


900 


900 


900 


900 


900 


Pressure 




900 






900 




























0. 20 






Porosity 


0. 35 


0. 35 


0. 25 


0. 25 


0. 25 


0. 25 


0. 25 


0. 25 


0. 25 


0. 25 


















0. 30 






Connate 


















0. 25 




Water 


0. 45 


0. 45 


0. 30 


0. 30 


0. 30 


0. 30 


0. 30 


0. 30 


0. 30 


0. 30 


Saturation 


















0. 35 




(m) 








0.1 




0.1 










Gas Cap Vol. 


0 . 1 


0. 1 


0. 5 


0. 5 


0. 1 


0. 5 


0 . 1 


0.1 


0.1 


0. 1 


Oil Vol. 








0. 9 




0. 9 










M 














0. 8M 








Viscosity 


1 . 0 


1 . 0 


1 . 0 


1 . 0 


1 . 0 


1.0 


1. 0M 


1 . 0 


1.0 


1 . 0 


Ratio 














1. 5M 








Case No. 


1 


1 


2 


2 


2 


2 


2 


2 


2 


2 


Figure 


4. 6 


5. 0 


6. 0 


7. 0 


8. 0 


9. 0 


10. 0 


11. 0 


12. 0 




No. 


4. 1 


5.1 


6.1 


7.1 


8.1 


9.1 


10. 1 


11. 1 


12.1 


none 










23. 5 




26. 0 


25 




22. 5 




% A N/N 


28. 6 


28. 6 


26. 5 


26. 5 


26. 0 


28. 0 


23. 6 


23. 6 


23. 6 


23. 6 










28. 6 




30.1 


21. 0 




25. 0 





V 
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produced directly from the gas cap. Since the computer program did not 
account for such gas production, an estimate of the production from the gas 
cap was made and production data were corrected accordingly. This non- 
ideal situation illustrates the flexibility of the program. Water influx is not 
included, but may be reflected by altering the relative permeability curve. 
Changing the input data requires some effort, but the effort is slight when 
compared with that required to carry out manually the Muskat Prediction. 

When the predicted results were analyzed and a set of parameters 
selected, it was necessary to consider whether or not the set of parameters 
chosen for Reservoir B was unique. 

Case I illustrated that no set of parameters, employed with porosity 
of 0. 30, 0. 35 and 0. 40, would satisfy observed production data. Consequently, 
when Case II variables were employed, it was quickly apparent that only a 
porosity of 0. 25 gave actual observed production in barrels of stock tank oil. 

N was not used. After the porosity was established, only the Kg/K 0 vs. S Q 
curve, which was positioned by S w 35 0. 30, gave a prediction which super- 
imposed on the observed production data. Since only this single combustion 
of S w and <t> resulted in observed production (% AN/N vs. Pressure) which 
superimposed on predicted production, the solution is a unique solution. These 
two variables were of primary importance because they were used to determine 
the observed % AN/N. The computed % AN/N is a function of (1 - S w ) times 
1/B 0 and does not reflect porosity at all. It follows that when an S w satisfies 
both the observed and predicted % AN/N, it is a unique value of S w . 
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CONCLUSION 

In this paper, the IBM- 701 digital computer solution of the Muskat 
equation was used to predict a solution gas drive reservoir. The conclusions 
drawn from this analysis are as follows: 

1. When using the Muskat Prediction to analyze a reservoir, one 
prediction is usually insufficient. 

2. A match between observed and predicted data is obtained only 
after several calculations. During these calculations, the 
different parameters were varied independently and in groups. 

3. Even when a match is obtained, the effect of variation of the 
remaining parameters must be; investigated to be reasonably 
sure that they will not destroy the match. 

4. The Pressure vs. % A N/N curve and the GOR vs. % A N/N 
curve could not be simultaneously matched with any set of 
reservoir parameters. The discrepancy between the observed 
and predicted GOR curves is attributed to gas production directly 
from the gas cap. 

5. The reservoir engineer's experience and knowledge of a reservoir 
are a necessary part of the procedure. He must decide when the 
parameters best specify the reservoir. Since only the ideal 
reservoir will conform exactly to the predicted solution, the 
reservoir engineer must decide when a reasonable approximation 
has been made and stop his analysis at this point. 
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NOMENCLATURE 



English 



Bo 

d 

f 

Gi 

Gp 



m 



n 

N p 

P 

Pab 
Pb 
P 

Q 
R 



con 



S g . 

S 0 i* 
S w 
U ■ 
X 
Y 



Greek 



r - 



a - 
- 

*- 



oil formation volume factor, reservoir barrels per stock tank barrel 
differential operator 

fraction of produced gas injected while producing below bubble-point 
pressure with declining pressure 

cumulative gas injected, thousands of cubic feet of gas measured at 
standard temperature and pressure per acre-ft. of formation 
cumulative gas produced, thousands of cubic feet of gas measured at 
standard temperature and pressure per acre-ft. of formation 
subscript, indicates initial or starting condition of pressure or 
saturation 

effective permeability to gas, md 
effective permeability to oil, md 

ratio of initial reservoir free gas volume to initial reservoir oil 
volume 

number of pressure or saturation steps 

cumulative oil produced, barrels of tank oil per acre-ft. of formation 

pressure, psia 

abandonment pressure, psia 

bubble-point (saturation) pressure, psia 

constant pressure, psia 

see Equation 3 

producing GOR, cubic feet of gas measured at standard temperature 
and pressure per stock tank barrel of oil 

solution GOR (gas solubility in oil), cubic feet of gas measured at 

standard temperature and pressure per stock tank barrel of oil 

gas saturation, fraction of pore space 

oil saturation, fraction of pore space 

oil saturation in the gas cap, fraction of pore space 

water saturation, fraction of pore space 

see Equation 3 

see Equation 3 

see Equation 3 



relative gas density, cubic feet of gas measured at standard 

temperature and pressure per cubic foot of gas measured at 

reservoir temperature and pressure 

difference ( P * PI - P 2 » S Q * S 2 - Sj ) 

gas viscosity, cp 

oil viscosity, cp 

porosity, fraction 
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CALCULATION OF INITIAL OIL IN PLACE 



Sand count of Figure 3: 2520 acre feet. 
Case I : 

N , 7758 (2520) (^ ) (1 - S w ) 

B 0 i 

1* B.R x1500, <t> • 0.35, S w * 0.45, B Q f 

N _ (7758) (2520) (0. 35) (0.55) 
1.185 

2. B. B* 900, i * 0. 35, S = 0. 45, B . « 

W Ol 

N = (7758) (2520) (0. 35) (0.55) 
1.113 

Case II: 

3. B. R* 1500, <6*0. 25, S w = 0. 30, B oi 

v . (7758) (2520) (0. 25) (0. 7) _ 

1.185 

4. B.R* 900, i * 0.25, S w = 0.30, B Qi 

KT _ (7758) (2520) (0. 25) (0. 7) 

IN = — * 3 

1.113 



« 1.185 

* 3. 18 x 10 6 ST bbls 
1.113 

= 3. 38 x 10 6 ST bbls 
= 1.185 

2. 88 x 10 6 ST bbls 

* 1.113 

.12 x 10 6 ST bbls 
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CALCULATION OF % A/V//V 



Pressure 

psia 


bbls of 
STO Produced 


Case I 

% 


Case II 

% 


GOR 

(actual) 


GOR 

(corrected) 


1100 


50, 000 


1. 48 


1. 6 


— 


— 


1000 


100, 000 


2. 96 


3. 2 


— 


— 


900 


175, 000 


5. 18 


5. 7 


600 


250 


800 


295, 000 


8. 75 


9. 5 


900 


500 


700 


410, 000 


12.1 


13.1 


1200 


800 


600 


490, 000 


14. 5 


15. 7 


1670 


1050 


Case I: N 


* 3. 38 million bbls of STO, <i> * 


• 35, 


= .45 




Case II: N 


■ 3. 12 million bbls. 


of STO, i 


= . 25, S 


w - .30 
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